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Soft elastic capsules which are driven through a viscous fluid undergo shape deformation coupled 
to their motion. We introduce an iterative solution scheme which couples hydrodynamic boundary 
integral methods and elastic shape eqnations to And the stationary axisymmetric shape and the 
velocity of an elastic capsule moving in a viscous fluid at low Reynolds numbers. We use this ap¬ 
proach to systematically study dynamical shape transitions of capsules with Hookean stretching and 
bending energies and spherical rest shape sedimenting under the influence of gravity or centrifugal 
forces. We And three types of possible axisymmetric stationary shapes for sedimenting capsules with 
fixed volume: a pseudospherical state, a pear-shaped state, and buckled shapes. Capsule shapes 
are controlled by two dimensionless parameters, the Foppl-von-Karman number characterizing the 
elastic properties and a Bond number characterizing the driving force. For increasing gravitational 
force the spherical shape transforms into a pear shape. For very large bending rigidity (very small 
Foppl-von-Karman number) this transition is discontinuous with shape hysteresis. The correspond¬ 
ing transition line terminates, however, in a critical point, such that the discontinuous transition is 
not present at typical Foppl-von-Karman numbers of synthetic capsules. In an additional bifurca¬ 
tion, buckled shapes occur upon increasing the gravitational force. This type of instability should be 
observable for generic synthetic capsules. All shape bifurcations can be resolved in the force-velocity 
relation of sedimenting capsules, where up to three capsule shapes with different velocities can occur 
for the same driving force. All three types of possible axisymmetric stationary shapes are stable 
with respect to rotation during sedimentation. Additionally, we study capsules pushed or pulled by 
a point force, where we always find capsule shapes to transform smoothly without bifurcations. 

PACS numbers: 47.57.ef,83.10.-y,46.70.De,46.32.-|-x 


I. INTRODUCTION 

On the microscale, the motion and deformation of 
closed soft elastic objects through a viscous fluid at low 
Reynolds numbers, either by a driving force or in a hy¬ 
drodynamic flow, represents an important problem with 
numerous applications, for example, for elastic microcap¬ 
sules [T], red blood cells nna, or vesicles moving in cap¬ 
illaries, deforming in shear flow [4], or sedimenting under 
gravity [S]. Another related system are droplets moving 
in a viscous fluid PIT]. 

The analytical description and the simulation of soft 
elastic objects moving in a fluid are challenging problems 
as the hydrodynamics of the fluid is coupled to the elastic 
deformation of the capsule or vesicle It is im¬ 

portant to recognize that this coupling is mutual: On the 
one hand, hydrodynamic forces deform a soft capsule, a 
vesicle, or a droplet. On the other hand, the deformed 
capsule, vesicle, or droplet changes the boundary condi¬ 
tions for the fluid flow. As a result of this interplay, the 
soft object deforms and takes on characteristic shapes; 
eventually there are dynamical transitions between dif¬ 
ferent shapes as a function of the driving force or flow ve¬ 
locity. Such shape changes have important consequences 
for applications or biological function. 
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In the following, we investigate the stationary shapes 
of sedimenting elastic capsules, which are moving in an 
otherwise quiescent incompressible fluid because of a ho¬ 
mogeneous body force, which can be either the gravita¬ 
tional or centrifugal force. We focus on capsules on the 
microscale with radii ~10 @m to 1 mm, a simple spherical 
rest shape, and sufficiently slow sedimenting velocities, 
such that the Reynolds numbers are small and we can 
use Stokes flow. 

We use the problem of sedimenting capsules to intro¬ 
duce a new iterative method to calculate efficiently ax¬ 
isymmetric stationary capsule shapes, which iterates be¬ 
tween a boundary integral method to solve the Stokes 
flow problem for given capsule shape [niHiisi and cap¬ 
sule shape equations to calculate the capsule shape in a 
given fluid velocity field. The method does not capture 
the dynamic evolution of the capsule shape but converges 
to its stationary shape. The method can be easily gener¬ 
alized to other types of driving forces apart from homo¬ 
geneous body forces. We show results for a driving point 
force but, in principle, the method can also be applied to 
self-propelled capsules. The methods allows us to inves¬ 
tigate instabilities and bifurcations of the final stationary 
capsule shape with high accuracy. 

Elastic capsules have a closed elastic membrane, which 
is a two-dimensional solid that can support in-plane shear 
stresses and inhomogeneous stretching stresses with re¬ 
spect to its rest shape and has a bending rigidity. The 
theoretical concept of an elastic capsule is rather general 
such that chemistry and nature provide many examples. 
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Synthetic capsules, for example, with polymer gel inter¬ 
faces, can be fabricated by various methods and with 
tunable mechanical properties and have numerous appli¬ 
cations in encapsulation and delivery [Bill]. There are 
also important biological examples of capsules, such as 
spherical (icosahedral) viruses, red blood cells, or artifi¬ 
cial cells consisting of a lipid bilayer and a cortex from 
filamentous proteins such as actin HMD, spectrin m, 
amyloid fibrils [m, or MreB filaments m- More gener¬ 
ally, the cortex of all animal cells consisting of the plasma 
membrane and the underlying actin (and spectrin) fil¬ 
ament network can be regarded as an elastic capsule. 
These types of capsules differ, however, in elastic prop¬ 
erties. 

Synthetic capsules from polymer gels typically have a 
spherical rest shape and sizes in the micrometer to mil¬ 
limeter range with a micrometer thick shell. Their shell 
materials can be described by isotropic elasticity of thin 
shells [nmn. Viruses are much smaller nanometer-sized 
objects but can also be described by thin shell isotropic 
elasticity [20]. Red blood cells acquire a nonspherical 
discoid rest shape under physiological conditions. Their 
elasticity is no longer isotropic but the bilayer mem¬ 
brane dominates the bending rigidity, whereas the spec¬ 
trin skeleton contributes a relatively small shear elastic¬ 
ity, while the bilayer membrane remains practically un- 
stretchable [mils]. For artificial or biological cells with 
an actin cortex the elastic properties are governed by the 
filamentous network and strongly vary with the mesh 
size: At large mesh size, the bending modulus is rela¬ 
tively large, resembling red blood cell elasticity, at small 
mesh size the bending modulus becomes small compared 
to the stretching modulus resembling the situation for a 
thin isotropic elastic shell. 

Sedimentation has already been studied for vesicles, 
both numerically [23U26] and experimentally |5|. As op¬ 
posed to elastic capsules, vesicles are bounded by a lipid 
bilayer, which is a two-dimensional fluid surface. For sed¬ 
imenting vesicles, pearlike and egglike shapes have been 
observed experimentally |5]. Numerically, depending on 
the initial configuration, pear shapes, banana shapes, or 
parachute shapes are found jlSl |21| . The banana shape 
exhibits circulating surface flows (33] [33]. At sufficiently 
high velocities, tethering instabilities occur both in ex¬ 
periments |S| and simulations |25j . 

Sedimentation has also been studied for red blood cells, 
which constitute a special type of soft elastic capsule, 
which is unstretchable and has a nonspherical biconcave 
rest shape. For sedimenting red blood cells shape tran¬ 
sitions also have been found. Early experiments on cen¬ 
trifuged red blood cells [371I3H] show shapes developing 
tails during centrifugation. In Ref. [391 a sequence of 
shapes from the biconcave rest shape to cup-shaped and 
bag-shaped cells has been observed as a function of sedi¬ 
menting time. Extensive multiparticle collision dynamics 
(MPCD) simulations on sedimenting red blood cell mod¬ 
els [30] produced a dynamic shape diagram, which ex¬ 
hibits tear drop shapes, parachute (or cup-shaped) blood 


cells, and fin-tailed shapes. 

Shape transitions of sedimenting spherical elastic cap¬ 
sules have not yet been studied although the spherical 
rest shape is prepared more easily in applications involv¬ 
ing synthetic microcapsules. Our iterative method allows 
us to completely characterize capsule shapes as a func¬ 
tion of their elastic properties and the driving force. The 
shapes depend on two dimensionless parameters: (i) the 
Foppl-von-Karman number describes the typical ratio of 
stretching to bending energy and characterizes the elastic 
properties of the capsule and (ii) the Bond number de¬ 
scribes the strength of the driving force relative to elastic 
deformation forces. We map out all stationary sediment¬ 
ing shapes in a shape dynamic diagram parameter plane 
spanned by Eoppl-von-Karman and Bond number. We 
can use this diagram to identify the accessible sediment¬ 
ing shapes for different types of capsule elasticity, such as 
isotropic shells, red blood cell elasticity, or semiflexible 
polymer network elasticity. 

Another important property of soft objects driven 
through a liquid is the relation between the driving force 
and the resulting velocity of the elastic objects. For 
strictly spherical sedimenting capsules this relation is 
the simple Stokes’ law. Shape transformations of a de¬ 
formable sedimenting object should reflect in qualitative 
changes or even bifurcations in the force-velocity curves. 
This important aspect has only been poorly studied for 
droplets, vesicles, and red blood cells. Some results 
have been obtained for quasispherical vesicles in Ref. [241 
MPCD simulations for red blood cells showed no specific 
signs of shape transitions in the force-velocity relations 
m, experimental measurements for vesicles [S] are also 
not precise enough to find such features. We will cal¬ 
culate the force-velocity relation for sedimenting elastic 
capsules with high accuracy in this paper, such that we 
can address this issue and identify shape bifurcations of 
sedimenting spherical capsules by their force-velocity re¬ 
lation. 


II. METHODS 

There are many simulation methods which have been 
applied to the dynamics of vesicles or capsules in viscous 
fluids, such as particle-based hydrodynamic simulation 
techniques such as MPCD [31], dissipative particle dy¬ 
namics (33), and Lattice Boltzmann simulations |33| or 
boundary integral methods [T|. These techniques simu¬ 
late the full time-dependent dynamics using a triangu¬ 
lated representation of vesicles or capsules; stationary 
states are obtained in the long-time limit. Here we intro¬ 
duce an iterative boundary integral method coupled to 
shape equations for axisymmetric shapes, which directly 
converges to stationary shapes without simulation of the 
real dynamics. Instead, we exploit axisymmetry to avoid 
triangulated representations and get an efficient iterative 
method based on boundary integrals for the fluid and 
shape equations for the capsule. 
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We limit ourselves to axisymmetric capsule shapes re¬ 
sulting also in axisymmetric viscous fluid flows. Then 
the elastic surface of the sedimenting capsule has to be 
at rest in its stationary shape without any surface flows. 
Therefore, also the fluid inside the capsule will be at rest 
(in the frame moving with the capsule). We will obtain 
the stationary capsule shape by solving a shape equa¬ 
tion in curvilinear cylindrical coordinates and solve for 
the Stokes flow of the surrounding fluid by applying a 
boundary integral method. In an iterative procedure we 
converge to the stationary shape, where the fluid forces 
onto the capsule and the boundary conditions for the 
fluid flow posed by the capsule shape are consistently in¬ 
corporated. We perform this iterative procedure for a 
given driving gravitational (or centrifugal) force and ob¬ 
tain the corresponding capsule sedimenting velocity from 
the condition of force balance with the total fluid force 
onto the capsule surface. 

A. Geometry 


n z 
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FIG. 1. (color online) Example of an axisymmetric shell with 
coordinates z, r, and ip. The shell is generated by revolution 
of the generatrix (orange thin line) that is calculated by the 
shape equations. 

We work in cylindrical coordinates to directly exploit 
the axisymmetry. The axis of symmetry is called z, the 
distance to this axis r, and the polar angle ip, see Fig. 
The shell is given by the generatrix [r(s),z(s)] which is 
parametrized in arc length s (starting at the lower apex 
with s = 0 and ending at the upper apex with s = L). 
The unit tangent vector to the generatrix at [r(s), ^(s)] 
defines an angle '0 via = (cos 0, sin 0), which can be 
used to quantify the orientation of a patch of the capsule 
relative to the axis of symmetry. 

B. Hydrodynamics 

We want to calculate the flow field of a viscous incom¬ 
pressible fluid around an axisymmetric capsule of given 
fixed shape. For the calculation of the flow field the de- 
formability of the capsule is not relevant, and the capsule 
can be viewed as a general immersed body of revolution 


B. For the calculation of the capsule shape and the de¬ 
termination of its sedimenting velocity, we only need to 
calculate the surface forces onto the capsule which are 
generated by the fluid flow. 

In the limit of small Reynolds numbers the stress ten¬ 
sor (T in a stationary fluid on which no external body 
forces are exerted is given by the stationary Stokes equa¬ 
tion [33] 

V • cr = 0. (1) 

In Cartesian coordinates, the stress tensor is given by 
o'ij = —pSij + ^pLCij, where p is the hydrodynamic pres¬ 
sure, /i the viscosity, and = (djUi + diUj) /2 is the 
rate of deformation tensor for a fluid velocity field u. 
Additionally, for an incompressible fluid, the continuity 
equation V • u = 0 holds. From the stress tensor one can 
infer the surface force density f = cr • n onto the capsule 
by means of the local surface normal n. 

In the rest frame of the sedimenting axisymmetric cap¬ 
sule and with a “no-slip” condition at the capsule surface 
dB, we are looking for axisymmetric solutions that have 
a given flow velocity at infinity and vanishing veloc¬ 
ity on the capsule boundary. In the laboratory frame, 
—u°° is the sedimenting velocity of the capsule in the 
stationary state. Therefore, has to be determined 
by balancing the total gravitational pulling force and the 
total hydrodynamic drag force on the capsule. 

In the laboratory frame we are looking for solutions 
with vanishing pressure and velocity at infinity. The 
Green’s function for these boundary conditions is the 
well-known Stokeslet (also called Oseen-Burgers tensor), 
that is, the fluid velocity u at y due to a point-force —Fp 
acting on the fluid at x, 

u(y) =-^G(y-x) • Fp (2) 

with the Stokeslet G whose elements are in Cartesian 
coordinates {x = |x|) 

= + ( 3 ) 

This is the basis of the boundary integral approach for a 
solution of the Stokes equation, where we want to find the 
correct distribution of point forces to match the “no-slip” 
condition on the capsule surface. 

Because of the axisymmetry we can integrate over the 
polar angle and find the velocity due to a ring of point 
forces with a local force density — f acting on the fluid, 
u(y) = ~8:^M(yiX) • f, with a matrix M, which can be 
calculated by integration of the Stokeslet G with respect 
to the polar angle. Switching from forces — f acting on 
the fluid to forces f acting on the capsule and integrating 
over all point forces on the surface, we find the general 
solution of the Stokes equation for an axisymmetric point 
force distribution on the surface of the body of revolution 

B, 

Wa(y) =-^^ds(x)M„,g(y,x)/^(x). (4) 
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Here Greek indices denote the component in cylindrical 
coordinates, i.e., a,/3 = r,z (u^ = fv = 0 for symmetry 
reasons). For these coordinates, the elements of the ma¬ 
trix kernel M are given in Appendix according to Ref. 
[HI The integration in Q runs along the path C given 
by the generatrix, i.e., the cross section of the boundary 
dB, with arc length s(x). This representation of a Stokes 
flow in terms of a single-layer potential (using only the 
Stokeslet and not the stresslet) is possible as long as there 
is no net flow through the surface of the capsule 

dA (u — Uq) • n = 0. (5) 

According to the “no-slip” condition we have (working 
in the laboratory frame) u = — u°° at every point y S dB 
on the surface. This results in the equation 

= ^ / ds(x)M„;3(y,x)/;3(x) (foryS^R), (6) 

which can be used to determine the surface force distri¬ 
bution. 

To numerically solve the integral for the surface force 
/(xi) at a given set of points {x^} [i = we 

employ a simple collocation method, i.e., we choose a 
discretized representation of the function fpijx) and ap¬ 
proximate the integral in ([^ by the rectangle method. 
This leads to a system of linear equations 

U = MF (7) 

with the “super-vectors” 

u = (u°°(yi)r, ^“(yi).^, • • • , M°°(yAr)r, U°^{yN)z) 

= ( 8 ) 

F = (/(xi)r, /(xi)^,. . . , /(xjv)r, /(xAr)^) (9) 

and a matrix (using numbered indices r=0, z=l for a 
more compact notation) 

M2i+a—l,2j-t-/3—1 — ('^(Xj-t-1) ■S (Xj ) ) Mq., 3 (y^, Xj ) . 

( 10 ) 

Due to singularities in the diagonal components of M, the 
two sets of points {x^}, {y^} must be distinct (or M other¬ 
wise regularized), so 2N points on the surface are needed, 
and the solution is the surface force at points {y^}. We 
validate our regularization in Appendix]^ For use in the 
shape equations, we then decompose this surface force 
into one component acting normal to the surface, the hy¬ 
drodynamic pressure p„, and another component acting 
tangential to the surface, the shear pressure Ps (see Fig. 

w 

We note that we usually restrict our computations to 
the bare minimum, i.e., the surface forces needed for the 
calculation of the capsule shape but, thereby, have all 
necessary information to reconstruct the whole velocity 
field in the surrounding liquid as shown in Fig. 


Throughout this paper, we use a “no-slip” condition, 
that is, the velocity directly at the surface of the im¬ 
mersed body vanishes in its resting frame. It is possible, 
however, to extend this boundary integral approach to in¬ 
corporate a prescribed velocity field on the surface in the 
resting frame of the capsule [8] . This will allow us to gen¬ 
eralize the approach to model active swimmers |351136] 
in future work. 


C. Equilibrium shape of capsule 

The sedimenting capsule is deformed by the hydrody¬ 
namic stresses from the surrounding fluid. We calculate 
the equilibrium shape of the capsule by a set of shape 
equations, which have been derived in Refs. [371 [351 We 
generalize these shape equations to include the additional 
fluid stresses on the capsule, which we obtain as the sur¬ 
face force density components fr and fz from the hy¬ 
drodynamic flow as described in the previous section. 
In order to make our approach self-consistent, the cap¬ 
sule shape B from which hydrodynamic surface forces 
are calculated has to be identical to the capsule shape 
that is obtained by integration of the shape equation un¬ 
der the influence of these hydrodynamic surface forces. 
This is achieved by an iterative procedure, which will be 
explained in the following section. 

The shape of a thin axisymmetric shell of thickness H 
can be derived from non-linear shell theory [391 SO] ■ A 
known reference shape [ro(so), Zo('So)] (a subscript zero 
refers to a quantity of the reference shape; sq G [0, Lq] 
is the arc length of the reference shape) is deformed by 
hydrodynamic forces exerted by the viscous flow. Each 
point [ro(so), 20 ( 50 )] is mapped onto a point [r(so)j 2;(so)] 
in the deformed configuration, which induces meridional 
and circumferential stretches. As = ds/dso and A,^ = 
r/ro, respectively. The arc length element ds of the de¬ 
formed configuration is = [r'(so)^ + z'{so)‘^]dsl. 

The shape of the deformed axisymmetric shell is given 
by the solution of a system of first-order differential equa¬ 
tions, henceforth referred to as the shape equations. Us¬ 
ing the notation of Refs. [371 [53] these can be written as 

r'(so) = As cosV' , z (so) = Xs simp , tp'{so)=XsKs 

r^(so) = As ^ cosip-i- Ks< 7 + Ps^ 

m(,(so) = As ^ costp - q 1^ 

q'{So) = As (^-KsTs - - ^COSIp p'^ . (11) 

The additional quantities appearing in these shape equa¬ 
tions are shown in Fig. and defined as follows: The an¬ 
gle Ip is the slope angle between the tangent plane to the 
deformed shape and the r axis, is the circumferential 
curvature, Ks the meridional curvature. The meridional 
and circumferential stresses are denoted by Ts and 
respectively; q is the transverse shear stress, p the total 
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FIG. 2. (Color online) Fluid velocity field around the four stationary capsule shapes under volume control, (a) pseudospherical, 
(b) pear shape, (c) strongly buckled, and (d) weakly buckled, in the capsule resting frame (left halves) and in the laboratory 
frame (right halves). The sedimenting force acts downwards, and arrows indicate the direction of flow. The arrows are 
equidistant in time, so a higher arrow density along a line represents a lower velocity. Additionally, the absolute value of the 
velocity is color coded in the background and the hydrodynamic surface forces are shown with arrows on the surface. The scale 
is identical in panels (a)-(d) and given by the scale bar of length Ro in (a). 


normal pressure, Ps the shear pressure, and I the external 
stress couple. 

The first three of the shape equations 0 follow from 
geometry, the last three equations express force and mo¬ 
ment equilibrium. For a derivation of these equations we 
refer to Refs. [37112S1 SSI In order to close the Eqs. ( [IT| ), 
additional geometrical relations 


^ ( 12 ) 

r as tq 


and constitutive stress-strain relations that depend on 
the elastic law governing the shell material are needed. 
The elastic law will be discussed in the following para¬ 
graph. 



FIG. 3. (Golor online) Left: Tensions and bending moments 
acting on a shell segment (thickness H). Right: Part of the 
generatrix (green thick line) with the decomposition of ex¬ 
ternal surface force densities into the coordinate components 
(/r-i /z) SI'S well as into normal and tangent components (pn, 
Ps). Also shown are the two principal curvature radii and 
K,~^ with corresponding osculating circles. 

The normal pressure 

P = Po + Pn - go^pz (13) 


(i.e., the pressure difference p = Pi^ — Pout between in¬ 
side and outside pressure), the shear-pressure Ps, and the 
stress couple I = PsH/2 (the fluid inside the capsule is 
at rest) are given externally by hydrodynamic and exter¬ 
nal driving forces; is the gravitational acceleration and 
Ap = Pin — Pout the density difference between the fluids 
inside and outside the capsule. Note that we measure the 
gravitational hydrostatic pressure —g^Apz relative to the 
lower apex, for which we can always choose 2 :( 0 ) = 0. 

The pressures Pn and Ps are the normal and tangential 
forces per area, which are generated by the surrounding 
fluid. The surface force density vector f = -|- frSr, 

which has been calculated in the previous section, can 
also be expressed in terms of p„ and ps, f = p„n — PsO^ 
(n and are normal and tangent unit vectors to the 
generatrix). Using the decompositions n = — cos + 
sin '06^ and eg = sin + cos with the slope angle 
Ip, we find: 


p„ = f • n = /r sin - /z cos ip (14) 

Ps =-i ■ Bs =-fr cosip - fzSinip. (15) 


This is also illustrated in Fig. 

All remaining quantities in the shape equations 0 
follow from constitutive (stress-strain) relations that de¬ 
pend on the elastic law, which will be discussed in the 
following paragraph. In addition to the elastic law, there 
might be global constraints. We consider here only one 
geometric constraint to the shape, namely a fixed vol¬ 
ume U = Vo (due to an incompressible fluid inside the 
closed capsule). We introduce the conjugated Lagrange 
parameter po, the static pressure difference between the 
interior and exterior fluids, which also enters the pressure 
p in Eq. (13). Additionally, the shape has to be in global 
force equilibrium, to which the velocity of the capsule is 
the conjugated parameter. 
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1. Elastic law and reference shapes 

Within this work we use a Hookean energy density 
with a spherical resting shape to model deformable cap¬ 
sules. As a measure for the deformation of the capsule 
we introduce the meridional and circumferential strains 

Cg — Ag 1 , Cip — \ip 1 (16) 

as well as the meridional and circumferential bending 
strains 


ATg = AgK, — K 


sO I 




(17) 


The elastic energy we use is Hookean, i.e., a quadratic 
form in these strains. More precisely, the surface en¬ 
ergy density w, which measures the elastic energy of an 
infinitesimal patch of the (deformed) surface divided by 
the area of this patch in the undeformed state, is given 
by 


Ws = \ 2i/ege,^ -f el) 

+ ^Eb {Kl + 2uKsK^ + Kl) 


(18) 


with the two-dimensional Young modulus Y 2 D, which de¬ 
fines the tension (energy per surface) scale, the bending 
modulus Eb and the (two-dimensional) Poisson ratio v 
(assuming equal Poisson ratios for bending and stretch¬ 
ing). The fundamental length scale is the radius i?o of 
the spherical rest shape (with Rq = KgjC^ and Lq = ttAq). 
Tension s an d bending moments derive from the surface 
energy ([l^ by Tg = X~^dw/des and ATg = X-^dw/dK^ 


(and two more analogous relations with indices s and 
ip interchanged), which gives the missing constitutive 
stress-strain relations for the shape equations (11). For 


Hookean elasticity and a spherical rest shape, the result¬ 
ing set of shape equations has been solved for a purely 
hydrostatic pressure p = po and pg = 0 in Ref. 133 

The Hookean elastic model for a spherical rest shape 
contains five parameters, Y 2 D, Aq, Eb, v, and El, char¬ 
acterizing capsule size and elastic properties. In Sec. 
|mg below, we will eliminate two parameters, Y 2 D and 
Ao) by choosing our scales of energy (or tension) and 
length appropriately. The two-dimensional Poisson ratio 
is bounded hy v ^ [—1; 1], and we reduce our parameter 
space by always using v = 1/2. This way our Hookean 
energy density has the same behavior for small stresses 
as the more complex Mooney-Rivlin functional [33] . The 
thickness H enters the shape equations only via the stress 
couple I and, thus, only weakly influences the resulting 
shape. Assuming that the shell can be treated as a thin 
shell made from an isotropic elastic material, the bending 
modulus is directly related to the shell thickness H by 


Thus, we have one remaining free parameter to change 
the elastic properties, the bending modulus Eb- 


The Hookean elastic energy law we use here is well 
suited to describe the deformation behavior of soft elas¬ 
tic capsules. For other systems different energies might 
be more adequate (e.g. Helfrich bending energy and a 
locally inextensible surface for vesicles). As pointed out 
above, a different choice of elastic law enters the formal¬ 
ism through the constitutive relations and, thus, does not 
require changes in our method on the conceptual level. 


2. Solution of the shape equations 


The boundary conditions for a shape that is closed and 
has no kinks at its poles are 

r(0) = r(Lo)=^(0)=7r-^(Lo)=0, (20) 

and we can always choose 2;(0) = 0. For the solution of 
the shape equations it is important that there is no net 
force on the capsule as the shape equations are derived 
from force equilibrium. The condition that the total hy¬ 
drodynamic drag force equals the total gravitational force 
determines the resulting sedimenting velocity. If hydro- 
dynamic drag and gravitational pull cancel each other in 
a stationary state, there is no remaining point force at 
the poles needed to ensure equilibrium and, thus. 


g(0) = qiLo) = 0. 


( 21 ) 


The shape equations have (removable) singularities at 
both poles; therefore, a numerical solution has to 
start at both poles requiring 12 boundary conditions 
(r, z, Ip, r, m, q at both poles), from which we know 7 [by 
Eqs. (20) and (21) and z(0) = Oj. The 5 remaining pa¬ 
rameters can be determined by a shooting method using 
that the solution starting at Sq = 0 and the one start¬ 
ing at So = Lq have to match continuously in the mid¬ 
dle, which gives 6 matching conditions {r, z,'ip,T,m,q). 
This gives an overdetermined nonlinear set of equations 
which we solve iteratively using linearizations. However, 
as in the static case m, the existence of a solution to 
the resulting system of linear equations (the matching 
conditions) is ensured by the existence of a first integral 
(see Appendix of the shape equation. In principle, 
this first integral could be used to cancel out the match¬ 
ing condition for one parameter (e.g., q), and, thus, we 
have 5 independent equations to determine 5 parame¬ 
ters, and the system is not genuinely overdetermined. 
We found the approach using an overdetermined system 
to be better numerically tractable, where we ultimately 
used a multiple shooting method including several match¬ 
ing points between the poles. Throughout this work we 
used a fourth-order Runge-Kutta scheme with step width 

Aso = 5 X IQ-^Ao. 

Using these boundary conditions, it is straightforward 
to see that the shape equations do not allow for a solution 
whose shape is the reference shape, unless there are no 
external loads {ps = p = I — 0). There will be solutions 
arbitrarily close to a sphere, which we call pseudospher- 
ical. 










7 


D. Iterative solution of shape and flow and 
determination of sedimenting velocity 


boundary condition 




external surface forces 


problem definition: 
reference shape 
gravity strength 
elastic constants 


FIG. 4. Iterative scheme for the solution to the problem of 
elastic capsules in Stokes flow as explained in the text. 


We find a joint solution to the shape equations and 
the Stokes equation by solving them separately and iter¬ 
atively, as illustrated in the scheme in Fig.|^ to converge 
to the desired solution: We assume a fixed axisymmetric 
shape and calculate the resulting hydrodynamic forces 
on the capsule for this shape. Then, we use the resulting 
hydrodynamic surface force density to calculate a new de¬ 
formed shape. Using this new shape we re-calculate the 
hydrodynamic surface forces and so on. We iterate until 
a fixed point is reached. At the fixed point, our approach 
is self-consistent, i.e., the capsule shape from which hy¬ 
drodynamic surface forces are calculated is identical to 
the capsule shape that is obtained by integration of the 
shape equation under the influence of exactly these hy¬ 
drodynamic surface forces. 

For each capsule shape during the iteration, we can 
determine its sedimenting velocity u = |u°°| by requiring 
that the total hydrodynamic drag force equals the total 
gravitational force, 


go^pVo = - [ ds(x)27tr(x)/2(x) (22) 

JC 

The integration runs along the path C given by the gen¬ 
eratrix, i.e., the cross section of the boundary dB, with 
arc length s(x). By changing u = |u°°| we can adjust 
the hydrodynamic drag forces on the right hand side to 
achieve equality for a given capsule shape during the iter¬ 
ation. During the iteration u will converge to the proper 
sedimenting velocity for the stationary state. 

The Stokes equation is linear in the velocity, and we 
can just scale the resulting surface forces accordingly, if 
we change the velocity parameter u°°. In this way, the 
global force balance can be treated the same way as other 
possible constraints like a fixed volume. Numerically, it 
is impossible to ensure the exact equality of the drag and 
the drive forces in Eq. (22). Demanding a very small 


residual force difference makes it difficult to find an ad¬ 
equate velocity, and a too-large force difference makes 


it impossible to find a solution with small errors at the 
matching points. Global force balance (22) is equivalent 
to the condition that the total force in axial direction 
vanishes, 


A(L) = —27T / ds r{p cos Ip + Ps sintp) = 0, (23) 

Jo 


see Eqs. (13) and (15) and Fig. Interestingly, the to¬ 
tal axial force is directly related to the existence of a 
first integral of the shape equations, which is discussed 
in the Appendixsee Eq. (C3). To monitor global force 
balance numerically, we use a criterion |A'(L)| < 10“^, 
which turns out to be a good compromise for a numerical 
force balance criterion. 


The iteration starts with a given (arbitrary) stress, e.g., 
one corresponding to the flow around the reference shape. 
For the resulting initial capsule shape, the Stokes flow is 
computed and the resulting stress is then used to start 
the iteration. If, during the iteration, the new and old 
stress differ strongly it might be difficult to find the new 
shooting parameters for the capsule shape and the right 
sedimenting velocity starting at their old values. To over¬ 
come this technical problem, we use a convex combina¬ 
tion (T = acTnew + (1 — a)o'oid of the two stresses and 
slowly increase the fraction a of the new stress until it 
reaches unity. The iteration continues until the change 
(monitored in pressure and velocity) within one iteration 
cycle is sufficiently small. In total, this allows for a joint 
solution of the shape equations and the Stokes equation 
to find the stationary shape and the sedimenting velocity 
at rather small numerical cost. 


If there are multiple stationary solutions at a given 
gravitational strength (see discussion of the shape dia¬ 
gram]^ in Sec. IIIB below), the iterative procedure ob¬ 
viously converges to only one shape. Which stationary 
shape is selected depends on the initial configuration. 
A stationary shape can be continuated to different pa¬ 
rameter values by slightly changing the control parame¬ 
ters and using the former stationary shape as new initial 
shape. 


In order to find all branches of stationary shapes we 
first generate capsule shapes for a fixed (artificial) flow 
field for the whole range of gravitational fields at a fixed 
bending rigidity. These serve as initial shapes from which 
we iterate until a stationary shape with the correct flow 
field is reached. At low bending rigidities, this procedure 
generates (all four) different classes of stationary shapes 
depending on the gravitational force. We then try to 
continuate all different classes of stationary shapes to the 
whole range of gravitational forces and, afterwards, to 
the whole range of bending rigidities. This allows us to 
obtain a full set of solutions and identify all different 
branches. 
















III. RESULTS 

In this section we present the results for stationary 
axisymmetric sedimenting shapes and stationary sedi¬ 
menting velocities as obtained by the fixed point itera¬ 
tion method. We focus on the sedimentation of capsules 
under volume control (V = Vq = const) in the main 
text and present additional results for pressure control 
(po = const) in Appendix Additionally, we show the 
results for a capsule that is driven (or pulled) by a local¬ 
ized point force rather than a homogeneous body force. 


(i) dimensionless bending modulus Eb or its inverse, the 
Foppl-von-Karman number 7 fvK = ^/Eb characterizes 
the elastic properties of the capsule (by the typical ratio 
of bending to stretching energy) and (ii) the Bond num¬ 
ber describes the strength of the driving force (relative to 
elastic deformation forces). If the capsule pressure is con¬ 
trolled rather than its volume, the dimensionless pressure 
Po provides a third parameter, see Appendix [E} Because 
stationary shapes are independent of time, they do not 
depend on the solvent viscosity. 


A. Control parameters and non-dimensionalization 


B. Shape diagram 


In order to identify the relevant control parameters and 
reduce the parameter space, in the remainder of the pa¬ 
per we introduce dimensionless quantities by measuring 
lengths in units of the radius Rq of the spherical rest 
shape of the capsule, energies in units of I 2 D.R 0 
tensions in units of I 2 D), and times in units of Ro/^/ 1 ^ 2 D- 
This results in the following set of control parameters for 
the capsule shape: 


1. Our elastic law is fully characterized by the dimen¬ 
sionless bending modulus or its inverse, the Foppl- 
von-Karman number [39) . 


Eb = 


Eb 

Y2dRI 


I 

7FvK 


(24) 


In general, elastic properties also depend on the 
Poisson ratio v but we limit ourselves to ;/ = 1/2 
as explained above. Using the thin-shell result 
(191 the dimensionless bending modulus Eb — 
(lI/9Ro)^ also determines the shell thickness H for 
an isotropic shell material. 


2. The sedimentation motion of the capsule is deter¬ 
mined by the strength of the gravitational (or cen¬ 
trifugal) pull (the gravitational force density), 
which we measure in units of I 2 D/Rq the follow¬ 
ing. The dimensionless gravitational force g then 
takes the form of a Bond number. 


The spherical rest shape has a dimensionless volume 
V = AnlZ. For sedimentation under volume control, we 
use this as the fixed volume Vq = 47t/3. In principle, 
it is also possible to prescribe volumes that differ from 
the rest shape volume. This can be done using pressure 
control, for which we present results in Appendix]^ For 
the static case, the only stationary shape at fixed volume 
is the strictly spherical shape, regardless of the bending 
modulus. This differs for a sedimenting capsule, which 
displays various shape transitions already at fixed vol¬ 
ume. 

We summarize our findings regarding the axisymmet¬ 
ric sedimenting shapes in the shape diagram Fig. in 
the {EB,g) plane of the two control parameters. We 
investigated stationary axisymmetric capsule shapes for 
gravity strengths (Bond numbers) up to g = Pmax = 2.5 
(usually using steps of Ag = 0.005) and bending moduli 
Eb = 10“^ to 10“^ corresponding to Foppl-von-Karman 
numbers 7 fvK = 10 to 1000. For high Bond and Foppl- 
von-Karman numbers, the iteration procedure becomes 
numerically more demanding (higher forces ask for finer 
discretization of the shape equations, smaller distances 
between solutions in parameter space ask for a more thor¬ 
ough fixed point search), but remains possible in princi¬ 
ple. We will discuss further below in detail, which pa¬ 
rameter regimes in the {EB,g) plane are accessible for 
different types of capsules and sedimenting driving forces. 

Overall, we identify three classes of axisymmetric 
shapes, which are shown in Fig. 


g = Bo = 


go^pRl 

F2D 


(25) 


where the characteristic elastic tension K 2 D is used 
instead of a liquid surface tension. 


3. The static pressure po within the capsule is mea¬ 
sured in units of ^ 20 /Ro in the following. 


1. Pseudospherical shapes, which remain convex and 
close to the rest shape. For pseudospherical shapes 
the velocity is still given by the result for a per¬ 
fect sphere, which is Stokes’ law u « F/6n « 2g/9 
(in dimensionless form), to a good approximation. 
Likewise, the pressure due to gravity is po « g{z) « 
9- 


4. The resulting sedimenting velocity u is measured in 
units of 12 d/m- 

We conclude that the resulting stationary sedimenting 
shapes of the capsule are fully determined by two dimen¬ 
sionless control parameters if the capsule volume is fixed: 


2. Pear shapes, which have a convex upper apex but 
develop an axisymmetric indentation at the side of 
the capsule. For high gravitational driving force 
(high sedimentation velocities) the pear shape re¬ 
sembles a tether with a high curvature at the upper 
apex. 
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FIG. 5. (Color online) Left: Stationary axisymmetric shapes for a sedimenting Hookean capsule with spherical reference shape 
in the parameter plane of the two control parameters, the dimensionless bending modulus Eb (the inverse of the Foppl-von- 
Karman number 7 fvk) and the dimensionless gravity g (the Bond number Bo). We find pseudospherical (blue), pear-shaped 
(red), and a pair of strongly (light green) and weakly (dark green) buckled shapes. The black line is a line of discontinuous 
shape transitions between pear and pseudospherical shape which terminates in a critical point beyond which a sphere transforms 
smoothly into the pear shape. The blue and red lines are spinodals indicating the limits of stability of the pseudospherical 
and pear-shaped shapes, respectively. Above the solid green line the pair of buckled shapes occur in a bifurcation; the dotted 
continuation of the solid line signals numerical difficulties in following this bifurcation line. To which solution the iteration 
procedure converges in regions with shape coexistence is controlled by the initial conditions. Right: Contours of all four types 
of shapes for Eb ~ 0.01, i.e., along the vertical black dotted line, in the range g = 0.7 to 1.7. The gravity strength g increases 
from left to right in steps of Ag = 0.2. 


3. Buckled shapes, which develop a concave axisym¬ 
metric indentation at the upper apex of the capsule. 
For fixed capsule volume, these buckled shapes al¬ 
ways occur in pairs of a weakly buckled shape with 
a shallow and narrow indentation and a strongly 
buckled shape with a deep and wide indentation at 
the upper apex. 

Our solution method allows us to identify all bifur¬ 
cations or transitions between these classes of shapes as 
shown in Fig.[^ Upon increasing the gravity g or decreas¬ 
ing the bending rigidity Eb, the stationary sedimenting 
pseudospherical shape transforms into a pear shape. The 
black line in Fig. is a line of discontinuous shape tran¬ 
sitions between pseudospherical and pear shape, which 
terminates in a critical point at 

Qc ~ 1.85 and Eb,c ^ 0.05. (26) 

At the discontinuous shape transition we find hysteresis 
with two spinodal lines (red and blue lines) indicating 
the limits of stability of pseudospherical and pear-shape 
shapes. For smaller bending rigidities Eb < Eb,c, the 
spherical shape transforms smoothly into the pear shape 
upon increasing gravity. We locate the transition lines 
between the sphere and pear shape by the condition of 
equal sedimenting velocity (for the discontinuous transi¬ 
tion) or maximal sedimenting velocity (for the smooth 
crossover) as explained in the next section. Discontinu¬ 
ous buckling transitions are also known from static spher¬ 
ical shells mi Eg, for example, as a function of the ex¬ 
ternal static pressure or an osmotic pressure [33] . That a 
line of discontinuous buckling transitions terminates in a 


critical point as a function of the elastic properties of the 
capsules is, however, an unknown phenomenon for static 
buckling transitions. 

Within the smooth crossover regime between sphere 
and pear shapes at small bending rigidities Eb < Eb^c, 
the pair of stationary buckled shapes occurs in a bifur¬ 
cation upon increasing the gravity strengths g above a 
critical threshold given by the solid green line in Fig. 

Also the spherical or pear shape remains a stable 
solution above the green line, such that we have three 
possible stationary sedimenting shapes in this parameter 
regime. Which of these shapes is dynamically selected in 
an actual experiment, depends on the initial conditions. 
The dotted continuation of the solid green line signals 
numerical difficulties in following this bifurcation line. 

The iterative method allows us to determine station¬ 
ary sedimenting shapes in the shape diagram Fig. very 
efficiently. In order to continuate the branch of pseudo¬ 
spherical and pear-shaped solutions for a fixed bending 
rigidity Eb = 0.01 in Bond number steps of Ag = 0.005 
in the range g = 0 , ..., 2 , we need a computing time of 
approximately 6 h on a four core Intel Xeon processor 
(3.7 GHz). 

C. Force-velocity relation 

All dynamic shape transitions reflect in bifurcations 
in the force-velocity relations. In Fig. we show the 
force-velocity relations between sedimenting velocity u 
and gravity g for all stationary axisymmetric shapes 
and for a color-coded range of bending rigidities Eb = 
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FIG. 6. (Color online) Force-velocity relations for sediment¬ 
ing capsules under volume control. The plot shows the veloc¬ 
ity relative to the (dimensionless) velocity mq = 2p/9 of a per¬ 
fect sphere with the same volume (Stokes’ law) as a function 
of the (dimensionless) total external force F = Angj?). The 
bending modulus Eb is color coded. Two distinct branches 
can be distinguished: The solution branches corresponding to 
pseudospherical and pear shapes start at zero force. A pair of 
solution branches for strongly or weakly buckled shapes starts 
at a nonzero threshold force. The branches for Eb = 0.01 are 
marked in black as a typical example. The capsule shapes 
1-6 in Fig. are realized at the points marked by circles. It 
is also apparent that the transition from pseudospherical to 
pear shapes is continuous for Eb ~ 0.05 and becomes discon¬ 
tinuous for Eb ~ 0.05. 


10“^ to 10“^. The total driving force F can be com¬ 
puted using the total drag from the fluid corresponding 
to the right hand side of Eq. (221 or as the product of the 
total mass difference and the sedimenting acceleration 
F = goApVo OT F = 47Tgf/3 (in dimensionless form) cor¬ 
responding to the left-hand side of Eq. (22). The leading 
contribution of the force-velocity relation is given by the 
force-velocity relation of a perfect sphere, which is Stokes’ 
law u = uq = F/6n = 2g/9 (in dimensionless form). To 
eliminate this leading contribution, we show the relative 
velocity u/uq in Fig.[^ By definition, the relative veloc¬ 
ity is also proportional to the sedimentation coefficient, 
which is defined as v/go and a standard quantity in cen¬ 
trifugation and sedimentation experiments |44) . 

Figure [^clearly shows two qualitatively different types 
of force-velocity curves. The force-velocity curves corre¬ 
sponding to pseudospherical or pear shapes with positive 
curvature at both apexes start at zero force, whereas the 
solution curves corresponding to the weakly and strongly 
buckled shapes bifurcate at a nonzero threshold force 
with a vertical tangent and have two branches: The 
weakly buckled shapes with a narrow indentation have 
a higher sedimenting velocity and correspond to the up¬ 
per branch, and the strongly buckled shapes have a lower 
sedimenting velocity and correspond to the lower branch. 

We note that, because of the g dependence of ug, the 
absolute sedimenting velocities u are increasing with in¬ 
creasing gravity for both buckled shapes, although u/ug 
is decreasing for the strongly buckled shapes. Above the 


threshold force for the buckled shapes, three different 
axisymmetric sedimenting capsule shapes can occur with 
different sedimenting velocities for the same gravitational 
driving force. Which of these shapes is dynamically se¬ 
lected in an actual experiment, depends on the initial 
conditions. 

The force-velocity curves for the pseudospherical or 
pear shapes allow us to detect how the transition be¬ 
tween spherical and pear shapes evolves into a discon¬ 
tinuous transition: Figure shows that these force- 
velocity curves develop a cusp close to the critical point 
for Eb ~ Eb,c — 0.05, whereas they remain smooth for 
Eb < Eb,c- For Eb > Eb,c, we find two overlapping and 
intersecting velocity branches, which corresponds to the 
characteristic hysteretic velocity switching in a discon¬ 
tinuous transition. In fact, we identify the discontinuous 
transition line in Fig. by these velocities: If the two 
branches coexist, then we localize the transition at the 
crossing of the two u{g) curves. For Eb < Eb,c, the ve¬ 
locity curve still exhibits a maximum, which we can use 
to define the crossover line between the pseudospherical 
and pear shape as it is shown in the shape diagram Fig. 

The discontinuous nature of the sphere-pear transition 
is confirmed by monitoring other quantities, such as the 
capsule area, as a function of the gravity g. As illustrated 
in Fig. the area clearly exhibits hysteretic switching at 
the transition. The pair of weakly and strongly buckled 
solutions appears in a bifurcation above a critical driving 
force. Also this bifurcation reflects in a corresponding bi¬ 
furcation of the capsule area as shown in the inset of Fig. 

m 



FIG. 7. (Color online) Reduced capsule area as a function 
of the dimensionless gravity strength g for a high bending 
modulus Eb = 0.08. There is an interval where pseudospher¬ 
ical and pear-shaped solutions with different area coexist with 
hysteretic switching upon increasing and decreasing g. Inset: 
For Eb = 0.01, the area changes continuously from pseudo¬ 
spherical to pear shapes. For the pair of weakly and strongly 
buckled solutions the area bifurcates at a critical g. 
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D. Transition mechanism 

Qualitatively, the shape transformations from a spher¬ 
ical shape into a pear shape or into buckled shapes are 
a result of buckling instabilities of a hydrodynamically 
stretched capsule at the upper part of the capsule. Upon 
increasing the gravitational force g the capsule acquires 
a higher sedimenting velocity and stretches along the z 
direction. This stretching is essential as it provides the 
excess area necessary for a shape transformation from a 
spherical rest shape, which has the minimal area for the 
fixed volume. 

The numerical results show that the fluid flow gen¬ 
erates a negative (compressive) contribution to the 
interior pressure at the bottom part of the capsule and 
a positive (dilatational) contribution at the upper part 
for all types of stationary shapes (see the orange arrows 
in Fig. [^. The gravitational hydrostatic pressure contri¬ 
bution —goApz, on the other hand, generates a negative 
contribution on the upper part as compared to a vanish¬ 
ing contribution on the lower apex (by choice of the z 
coordinate). Together with the positive volume pressure 
Pq, this results in a compressive negative total normal 
pressure at the upper part and dilatational positive total 
pressure at the lower part. Both for pear and buckled 
shapes, the negative total normal pressure at the upper 
part of the capsule, which is mainly caused by the grav¬ 
itational hydrostatic pressure contribution, is the reason 
to develop indentations. In the transition from pseudo- 
spherical to pear-shaped capsules an axisymmetric inden¬ 
tation develops at the side of the capsule (see the shapes 
in the middle row in Fig. [^, because for an elongated 
stretched capsule, this part of the capsule has lower cur¬ 
vature and, thus, less stability with respect to buckling. 
At higher driving force g and sedimenting velocities, the 
indentation can also form at the upper apex (see top and 
bottom rows in Fig. [^, and the pair of buckled configu¬ 
rations become stable solutions. 

An approximative limit for the stability of the spheri¬ 
cal shape can be given in terms of the classical buckling 
pressure [39l|4Tl|42] Pc = —4\/ Eb in dimensionless units. 
We first note that the height of a sphere is z„iax ~ 2 
and that the static pressure contribution needed to main¬ 
tain a constant volume against the hydrostatic pressure 
is po « gz, where z « 1 is the z-component of the cen¬ 
ter of mass. This gives an effective compressive pressure 
Po — gzmax ~ —g at the upper apex (the compressive hy¬ 
drodynamic contribution > 0 is smaller and can be 
neglected compared to relative to hydrostatic and static 
pressure). If this compressive pressure exceeds the clas¬ 
sical buckling pressure, i.e., g > \pc\ = 4-\/ Eg, then the 
spherical shape is unstable with respect to indentations 
at the upper side of the capsule, which leads to the pear 
shape. A parameter dependence 



describes well the boundary between spherical and pear 


shape in the shape diagram Fig.[^ 

The termination of the line of discontinuous transitions 
between sphere and pear shapes terminating in a critical 
point is a result of the deformation of the capsule by the 
fluid flow prior to the shape transition: Increasing the 
gravity g gives rise to a hydrodynamic stretching of the 
upper part of the capsule. The smaller the dimensionless 
bending rigidity As, the smaller is the meridional curva¬ 
ture Ks in the upper part. For soft capsules, the merid¬ 
ional curvature vanishes before the effective compressive 
pressure becomes comparable to the buckling pressure, 
(see, for example, the second shape in the middle row in 
Fig. I for Eb = 0.01). For a flat shell segment, how¬ 
ever, buckling becomes continuous and does not require 
a threshold normal pressure. Rigid capsules, on the other 
hand, remain curved upon increasing the gravity g up to 
the discontinuous buckling induced by the normal pres¬ 
sure in the upper part of the capsule. 

For the bifurcation line of the two buckled shapes in 

the shape diagram Fig. [5] we find an approximately lin- 
-1,7 

ear dependence g oc \/ Eb from fitting our numerical 
results. Currently, we have no simple explanation for this 
result. 


E. Rotational stability 

We always assumed axisymmetric shapes that perform 
a rotation-free sedimentation, i.e., stability with respect 
to rotation or tilt of the shape. The question of whether 
the sedimentation motion of an axisymmetric rigid body 
is rotationally stable can be reduced [34] to the question 
whether the so-called center of hydrodynamic stress lies 
above or below the center of mass; both points have to 
be on the symmetry axis. In general, sedimenting bod¬ 
ies will tilt in a way that aligns the connection vector 
of the centers of mass and hydrodynamic stress with the 
direction of gravity [IS] or start to rotate. The center of 
hydrodynamic stress is the point about which the trans¬ 
lational and rotational motions decouple. 

In Appendix [^ we present a stability argument that 
only uses the information available to us. The result of 
this analysis is exemplarily shown in Fig. Ufor Eb = 0.01. 
We find that the center of hydrodynamic stress zq lies 
almost always below the center of gravity Zc.m. for all 
three classes of shapes, pseudospherical, pear-shaped, 
and buckled such that the shapes are linearly stable 
against out-of-axis rotations. Only solutions in the pear- 
shape class for very high gravity forces g > 2.2 are unsta¬ 
ble. These shapes exhibit a very pronounced tether-like 
extrusion. 


F. Accessible parameter space 

With the shape diagram in Fig. [^we have a complete 
overview of theoretically possible shapes and transitions 
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FIG. 8. (Color online) Linear analysis concerning rotational 
stability for Eb = 0.01. We show the (dimensionless) dif¬ 
ference zo — 2c.m. between the center of hydrodynamic stress 
2o and the center of gravity 2c.m. as a function of the di¬ 
mensionless gravity strength g for all three solution branches 
[(blue) pseudospherical, (red) pear-shaped, (green) buckled]. 
For 2 o — 2c,m. < 0, i.e., below the gray line, the shape is lin¬ 
early stable against out-of-axis rotations. Only solutions with 
a very pronounced tether extrusion at high g are unstable. 


between them in the plane of the two control parameters, 
the dimensionless bending rigidity Eb or its inverse, the 
Foppl-von-Karman number 7 fvK = ^IEb, and the Bond 
number Bo = g. There are three limitations, however, to 
the accessible parameter range. 

First, our analysis is limited to the regime of low 
Reynolds numbers because we used the Stokes equation 
to describe the fluid flow. The Reynolds number for the 
sedimenting capsule is given by Re = vpRo/g,. To lead¬ 
ing order the sedimenting velocity is given by Stokes’ law 
V = goApV/dnpRo resulting in Re = 2goAppiig/9^^. A 
criterion Re < 1 limits the accessible Bond numbers to 
9 < 5max, Re ~ 9M^/2p.Ro^2D ■ For water as solvent we 
have gmax, Re 5 X Rq^Y^ N. This can be easily 
increased by a factor of 100 in a more viscous solvent. 
For a capsule of size i?o ~ lO/rm in water, the criterion 
g < 5max, Re corresponds to a condition that the accelera¬ 
tion should remain smaller than 5000 times the standard 
gravity. 

Second, the sedimenting force is very small for gravi¬ 
tation; it is much larger for a centrifuge but also modern 
ultracentrifuges are limited to accelerations about 10® 
times the standard gravity [44]. With a typical density 
difference of Ap = 10^ kg/m® (10% of the density of wa¬ 
ter) we find that the Bond number is limited in available 
centrifuges to g 7 Pmax, cf ^ IQ^RlY-^ N/m®. 

Both of these constraints give an upper bound for the 
accessible Bond number range. The constraints are size 
dependent, however: For small capsule sizes i?o ^ 1 
for example, for viruses, the centrifuge constraint is more 
restrictive, whereas for larger capsule sizes Ro // 1 
the low-Reynolds-number constraint is more restrictive. 
The latter case is typical for all synthetic capsules, red 
blood cells, or artificial cells consisting of a lipid bilayer 


and a filamentous cortex, which all have capsule sizes in 
the range i?o 10 to 1000pm. 

Third, it is not possible to realize arbitrarily large di¬ 
mensionless bending rigidities Eb = Eb/Y 2 uRo or small 
Foppl-von-Karman numbers in experimentally available 
capsules. For isotropic elastic shell materials the rel¬ 
ative shell thickness determines this parameter, Eb ^ 
(iL/i?o)^) see Eq. (191, and values Eb ^ 0.01 are difficult 
to realize with synthetic capsules. Viruses also follow this 
law (with a shell thickness H ~ 2nm), resulting in Foppl- 
von-Karman numbers 7 fvK ^ 60 or Eb ^ 0.02 with Eb 
decreasing with increasing virus size |20j : small viruses 
with Rq ^ 15nm realize the largest values of Eb- 

Red blood cells also have similar values ot Eb al¬ 
though they are elastically strongly anisotropic because 
they consist of a lipid bilayer and a spectrin cytoskeleton. 
The liquid lipid bilayer dominates the bending modulus 
Eb SO/csT, whereas the shear modulus is determined 
by the filamentous spectrin cortex or skeleton. Filamen¬ 
tous networks typically have low shear moduli, for exam¬ 
ple G 2 D ~ 10“®N/m for the spectrin skeleton of a red 
blood cell Uni- Using V 2 D 4G2d for the red blood cell, 
because the lipid bilayer has a much higher area expan¬ 
sion modules K 2 B ~ 0.3N/m, we obtain Eb 0.005 for 
a size Rq ~ 3pm. Larger values oi Eb could be realized 
in artificial cells consisting of a lipid bilayer and a cor¬ 
tex from filamentous proteins for low filament densities 
in the cortex. These estimates of Eb for red blood cells 
or artificial cells might, however, be misleading because 
the relevant elastic modulus in the elastic energy is 
actually V 2 d /(1 — rather than F 2 D, which is used in 
the standard definition of the Fbppl-von-Karman number 
and Eb, see Eq. (24). Because of their lipid bilayer, red 
blood cells or artificial cells are in the unstretchable limit 
K 2 B ^ G 2 D) where v approaches unity, and we obtain 
V2 d/(1 — 1 ^^) ~ K 2 d- Using this modulus in a definition 
of Eb (i.e., in the non-dimensionalization) leads to much 
smaller values Eb ~ 10“®. 

The discontinuous sphere-pear transition happens be¬ 
yond a critical point, Eb > Eb,c — 0.05, see Fig. ^ 
Therefore, this peculiar transition is not accessible for 
typical synthetic capsules. To observe this transition an 
elastically very anisotropic capsule shell material with 
high bending and low stretching moduli would be needed. 
The transition into buckled shapes, however, should be 
well accessible with generic soft synthetic capsules. 

The softness of the shell material is also important if 
high Bond numbers Bo ^ 1 are to be realized in the 
shape diagram Fig. In the presence of the above 
low-Reynolds-number constraint, very soft shell materi¬ 
als with a small two-dimensional Young modulus Y 2 D 
5 X lO^'^N/m for a capsule of size i?o ~ 10pm are needed 
to reach such Bond numbers in water. Typical soft syn¬ 
thetic capsules such as OTS capsules have much larger 
moduli Y 2 B O.lN/m (38). Red blood cells or artificial 
cells with shear moduli governed by a soft filamentous 
network have small moduli Y 2 D ~ 10“®N/m, which will 
allow us to reach high Bond numbers. Again, these es- 
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timates might be misleading because the relevant elastic 
modulus in the elastic energy is actually F 2 d /(1 — 
rather than 120 - If this modulus is used in the nondi- 
mensionalization and the definition of the Bond number 
in Eq. (25), Bond numbers become much smaller because 
red blood cells or artihcial cells are hardly stretchable. 


IV. LOCALIZED DRIVING FORCES 

There are other possible external driving forces and 
self-propulsion mechanisms [Ml HZ]- One extreme case 
is a very localized external force. We study this case 
by employing a driving pressure that acts on a small 
patch on one apex (we used the criterion sq ^ 0-1 or 
L — sq < 0.1, respectively; even more localized pressures 
would have to be higher and, thus, require smaller inte¬ 
gration steps). For a pushing force, i.e., a force that acts 
on the stern of the capsule, this gives rise to strongly in¬ 
dented shapes. They have a higher drag than a spherical 
shape such that the driving force has to be higher than in 
Stokes’ law u « FjQn (in dimensionless form) resulting 
in Fdrive > 6™, whereas for pulling forces (acting on the 
bow) the drag is lower as compared to a sphere. 

We show force-velocity relations with some exemplary 
shapes for these two types of external driving forces in 
Fig.[9l As can be seen clearly seen in the force-velocity 
relations, there are no shape bifurcations for elastic cap¬ 
sules driven by localized forces. 

One way to apply a localized force experimentally is 
to attach beads to the capsule that can be manipulated 
via optical tweezers, as it has been done (for pulling and 
pushing types of forces) with giant vesicles in Ref. |48l 


V. DISCUSSION AND CONCLUSION 


We introduced a new iterative solution scheme to find 
stationary axisymmetric shapes and velocity of a de¬ 
formable elastic capsule driven through a viscous fluid by 
an external force at low Reynolds numbers. We focused 
on homogeneous body forces, i.e., sedimenting capsules 
in gravitation or in a centrifuge, but also demonstrated 
that other force distributions, such as localized forces, 
can be studied. 

The iterative method is sufficiently accurate and fast 
to identify all branches of different shapes and to resolve 
dynamic shape transitions. We find a rich dynamic bifur¬ 
cation behavior for sedimenting elastic capsules, even at 
fixed volume, see Fig. with three types of possible ax¬ 
isymmetric stationary shapes: a pseudospherical state, a 
pear-shaped state, and buckled shapes. The capsule can 
undergo shape transformations as a function of two di¬ 
mensionless parameters. The elastic properties are char¬ 
acterized by a dimensionless bending rigidity Eb [the in¬ 
verse Foppl-von-Karman number, see Eq. (24)], and the 
driving force is characterized by the dimensionless Bond 
number [see Eq. (25)]. 



FIG. 9. (Color online) Force-velocity relations for a cap¬ 
sule with fixed volume driven by (a) a localized pulling force 
and (b) a localized pushing force. The plots show the velocity 
relative to the (dimensionless) velocity uq = F/Gn of a hard 
sphere with same volume as a function of the (dimension¬ 
less) total external force. The insets show exemplary capsule 
shapes for Fb = 0.01. 


The transition between pseudospherical and pear 
shape is a discontinuous transition with shape coexis¬ 
tence and hysteresis for large bending rigidity or large 
Bond number. The corresponding transition line termi¬ 
nates, however, in a critical point if the bending rigidity 
is lowered, a phenomenon which is unknown from dis¬ 
continuous static buckling transitions. Parameter esti¬ 
mates show that this transition cannot be observed with 
standard synthetic capsules because it requires a rather 
large bending rigidity or low Foppl-von-Karman number, 
which cannot be realized for thin {H/Rq < 0.1) shells of 
an isotropic elastic material. Observation of this discon¬ 
tinuous transition requires a material, which combines 
high bending rigidity while it remains highly stretchable. 

We find an additional bifurcation into a pair of buck¬ 
led shapes at small bending rigidities upon increasing 
the gravitational force. This bifurcation should be ex¬ 
perimentally accessible with synthetic capsules in a cen¬ 
trifuge. The bifurcations are caused by hydrodynamic 
stretching, which increases a compressive hydrostatic 
pressure at the upper apex and leads to buckling-like 
transitions. 

All shape bifurcations can be resolved in the force- 
velocity relation of sedimenting capsules, where up to 
three capsule shapes with different velocities can occur 
for the same gravitational driving force, see Fig. In 
an experiment, these shapes are selected depending on 
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the initial conditions. We also confirmed the stability of 
all axisymmetric shapes with respect to rotation around 
an axis perpendicular to the axis of symmetry, see Fig. 
1^ We find more stationary buckled shapes if we con¬ 
sider sedimenting capsules where we vary the volume by 
changing the pressure as it is explained in more detail 
in Appendix]^ These additional shapes are higher order 
buckled shapes with more indentations. 


It is instructive to discuss our results in compar¬ 
ison with sedimenting vesicles [53H2S] and sediment¬ 
ing red blood cells m- Vesicles are bounded by a 
two-dimensional fluid lipid bilayer rather than a two- 
dimensional solid membrane or shell and are virtually 
unstretchable, i.e., they have a fixed area. Therefore, 
shape transitions from a spherical shape with the mini¬ 
mal area for given volume require excess area. For the 
sedimenting capsules this additional area is generated by 
hydrodynamic stretching. For vesicles, excess area can 
be “hidden” in thermal fluctuations. A spherical vesicle 
without such excess area cannot undergo any shape tran¬ 
sitions during sedimentation. The pear shape for elastic 
capsules at high Bond numbers is similar to the “tether” 
formation by sedimenting vesicles [25] . However, there 
are two important differences. First, the stretching en¬ 
ergy of capsules penalizes deformations from the rest¬ 
ing shape (whereas vesicles have a liquid surface that 
allows any shape with the correct area), which leads to 
less extreme and pronounced tethers. Moreover, cylin¬ 
drical tethers are an actual equilibrium shape of a mem¬ 
brane under tension, which can coexist with a spherical 
vesicle [49]. This is not the case for an elastic capsule. 
Secondly, we do not find a bulge or “droplet” forming 
at the upper end, by which we mean an increase in the 
width of the extrusion near the end of the tether. This 
is rooted in the different bending energies: Using a Hel- 
frich bending energy that is quadratic in the mean cur¬ 
vature wb (ks + the negative meridional curva¬ 

ture {ks < 0) needed to form a drop at the end of the 
tether is energetically favorable, whereas it is not with 
the Hookean bending energy (18) we employ. Moreover, 
elastic capsules cannot develop shapes with circulating 
surface flows because of their solid membrane. There¬ 
fore, there is no analog of the banana shape with surface 
flows that has been found in Refs. HSlIMlfor vesicles. 


Also sedimenting red blood cells exhibits different 
shapes; in MPCD simulations transitions between tear 
drop shapes, parachute (or cup-shaped) blood cells and 
fin-tailed shapes [30] have been found. Red blood cells 
have a non-spherical discocyte rest shape, which already 
provides some excess area as compared to the minimal 
area of a sphere. Therefore, stretching is not essential for 
dynamic shape transitions. The initial discocyte shape of 
red blood cells also gives rise to a strong tendency to tilt, 
which is absent for our initially spherical capsules. The 
parachute shape also tilts by almost 90° during transfor¬ 
mation into a tear drop [30] . Interestingly, all transitions 
observed in Ref. [sni did not show any signatures in the 
force-velocity relation or the bending or shear energies. 


For spherical capsules we find clear signatures of all tran¬ 
sitions in the force-velocity relation in Fig. [^ It remains 
to be clarified in future work whether the MPCD tech¬ 
nique is not accurate enough or exhibits too many fluctu¬ 
ations to resolve these signatures or whether the dynamic 
transitions of sedimenting red blood cells qualitatively 
differ, for example, because tilt plays such a prominent 
role. 

Finally, we studied capsules driven by localized forces, 
such as a pushing or pulling point-like force. For such 
localized driving forces, we always find shapes to trans¬ 
form smoothly without any shape bifurcations, see Fig. 


We plan to extend our iterative method to actively 
swimming capsules in future work |36j . One possi¬ 
ble swimming mechanism, for example in squirmer type 
swimmer models |50j , is a finite velocity field at the cap¬ 
sule surface in its resting frame. Our solution technique 
using single-layer potentials remains still applicable as 
long as there is no net velocity normal to the surface 
(which indicates a fluid source or sink within the capsule 
and is, thus, not relevant). 
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Appendix A: Single-layer potential solution of the 
Stokes equation in an axisymmetric domain 


To make the paper self-contained, we give the analyt¬ 
ical expressions for the elements of the matrix kernel M 
that relates the forces and the velocities on the surface 
of the capsule via 


= ds(x)M„; 3 (y,x)/; 3 (x) (foryeaU), 

^ .(Al) 

see Eq. Q. The derivation of this matrix kernel is briefly 
outlined in the main text and can also be found in Ref. 

IHl 

As we are working in cylindrical coordinates (r,z,ip) 
and the problem is axisymmetric, the matrix has four 
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elements (a, l3 = r, z) which can be explicitly written as 
M,,(y,x) = 2fc + (A2) 


Mzr(y,x) = 


kz 


z 

Vv 


{F-{rl-rl+z^)E) 


(A3) 


M,,(y,x) = -k-J- {F + {rl - - z^)E) (A4) 

Ty y Ty 


^(y>x) = 


' x' y y 


^-F[^rl + rl + 2z^)F 


-{2z^ + iz^irl + rl) + (r^ - rl)^)E] (A5) 

using the abbreviations 

(A6) 
(A7) 


Z = Zx — z. 


X '^y 


k^ = 
E = 




£2 + (r^ + ry)2 

E{k) 

F = K{k) 


(AS) 

(A9) 


and the complete elliptic integrals of the first and second 
kinds m, 



dc 


K{k) = 

Jo Vl — k? C0s2 X 

(AlO) 


r'^/2 _ 


E{k) = 

/ dc V 1 — cos^ x 

Jo 

(All) 


For the numerical computation of M elliptic integrals of 
the first and second kind we use polynomial approxi- 
mants as given in Ref. IMl 


Appendix B: Validation of numerical method 

We validate our algorithm and, in particular, our nu¬ 
merical treatment of singularities occurring in the bound¬ 
ary integral approach in Eq. Q by comparison with exact 
analytical results for the Perrin factors for the total drag 
of a spheroid with semiaxes a in radial and b in axial 
direction [34l [52]. A spheroid is an ellipsoid with two 
degenerate semiaxes a in radial direction, i.e., it can be 
parametrized by 


( a cos 0 sin (^ \ 

asin0sin(/? (Bl) 

b cos 9 j 

with the polar angle 0 < (/? < 27T and the parametric 
latitude 0 < 0 < tt. 

The Perrin factor S is the ratio of the total drag force 
onto a spheroid with semiaxes ratio ^ = b/a (for transla¬ 
tion along the axial a direction) and the drag of a sphere 
of the same volume moving with the same velocity. S is 



FIG. 10. (Color online) Perrin factor H, i.e., drag force onto 
a spheroid in a viscous fluid relative to that onto a sphere 
of same volume moving with the same velocity as a function 
of the ratio of the spheroid semiaxes ^ = b/a. The theoreti¬ 
cal result (B2 l due to Perrin is shown as a solid (black) line 
and our numerical data as (green) crosses. The relative dif¬ 
ference is less than 1/1000 for all values, which validates our 
numerical approach. 


known analytically (from calculating the stream function 
in ellipsoidal coordinates) [52] . 


3 

4 


1 


g _ 

VSp-l(-5p-l-(5p-l-l) artanh({p ^)) 

ri 




for ^ > 1 

(B2) 

for ^ < 1 


with 


^ ' (B3) 

Co = (B4) 


for both prolate (^ > 1) and oblate (^ < 1) spheroids. 
We compare our numerical results to this analytical ex¬ 
pression in Fig. [10|and find excellent agreement. 


Appendix C: First Integral of the shape equations 


Inspired from the known first integral in the static 
problem (cf. Eqs. (17) and (22) in Ref. [ST]) we make the 
following ansatz for a first integral of the shape equations: 

U(s) = 27tr (q cos i/) + TsSinil)) + X = const. (Cl) 

We are also assuming that the pressure p and the shear 
pressure Ps can be written as functions solely dependent 
on the arc length s. The calculation is rather straight¬ 
forward: We differentiate using the shape equations 0 
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and the geometric relations (12) and get 
0 = U'{s) 

= 2nq cos^ ip — 2nrKsq sin ip 
+ 2T[r cosip (^{ — KgTs — KipTtp — cosip- + 

+ 2nTs cos Ip sin ip + 2nrKsTs cos ip 

+ 2nrsinip ^osip——— + Kgq + Ps^ + X' 

= 2nrp cos ip + 2nrps sin ip + X' = Q (C2) 


In the second-to-last step most terms cancel each other 
out. Thus, we find an ordinary differential equation for 
X, which we can integrate directly: 

X' = —2nr{p cos ip + Ps sin ip) 

X = —2n dx r{p cos Ip + Ps sin Ip) (C3) 
Jo 

Inspecting the behavior at s = 0 we then deduce U(s) = 
U(0) = 0 or 

0 = 27xrq cos Ip + 2nrTs sin ip + X. (C4) 


The first integral U{s) = 0 is indeed a generalization of 
the static quantity, as we see by setting ps = 0 , p = const, 
which gives 


X = —2np / dxr cosip = —npr^, (C5) 
Jo 


and, thus, Eq. (C4| becomes equivalent to Eq. (17) in 
Ref. [371 

The first integral is related to the global axial force 
balance of the capsule, which can be seen by considering 
s = L, where the U{s) = U{L) = 0 gives X{L) = 0. The 
quantity X contains the contribution to the net force in 
the z direction and thus a shape with the desired features 
(namely g = 0 at the apexes) must have X{L) — 0 and, 
thus, be in global force balance in axial direction. 


Appendix D: Rotational stability 

We want to study the change in the torque due to an 
infinitesimally small rotation R of the velocity and the 
gravitational force vector. Geometrically, this is equiv¬ 
alent to rotating the capsule. A change of the velocity 
boundary condition Uq —>■ Uq = Ruq leads to new hydro- 
dynamic surface forces £—>■£'. As we did not change the 
capsule we can employ the reciprocal theorem 

J dAu-f' = J dAug-f (Dl) 

and find f' = (R^)“^f = Rf. The new hydrodynamic 
torque is, thus, 

Tu' = y dA (r X (Rf)) = RjdA ((R^^r) x f) . (D2) 


We limit ourselves to a linear analysis, that is, we perform 
an infinitesimal rotation around the y axis by an angle 
da, i.e., R^^ = I ± daJy. Here I is the unit matrix (/^ = 
5ij) and Jy generates a rotation around the y axis with 
~ ~ Then, to linear order in da, 

Gy ■ Tu' = -27tda J dsr Q^/r + {z - zo)/z^ ■ (D3) 

Here zq gives the pivot point of the rotation. Since we 
know that the centers of mass and hydrodynamic stress 
have to lie on the symmetry axis, studying rotations 
about points on this axis suffices. Likewise, we can cal¬ 
culate the gravitational contribution to the new torque, 
which is (again to linear order in da) 


By • Tg'= Ttdag J dsr‘^{z — zo). (D4) 

Equating these two gives the pivot point of the 
marginally stable infinitesimal rotation, i.e., the center 
of hydrodynamic stress. 

This argument does not account for a change of the 
capsule shape due to the altered stresses upon rotation, 
which means that a seemingly unstable capsule might be 
stabilized axially by a “wobbling” deformation leading to 
different hydrodynamic torques. In principle, this might 
also work the other way round (as for semiflexible cylin¬ 
drical rods [53]), but we think it is less relevant here, 
because of the typically large radius of curvature at the 
lower apex. 


Appendix E: Pressure control 


For volume control, we fixed Vq = 47 t/3 to the volume 
of the spherical rest shape and only explored the two pa¬ 
rameters Eb and g. The pressure po served as Lagrange 
parameter to achieve the fixed volume Vq. To get a better 
understanding of the variety of shapes that are in princi¬ 
ple possible, we study capsules with volumes other than 
that of the reference shape by controlling and changing 
the Lagrange parameter pq. Thus, the explorable param¬ 
eter space now consists of three parameters Eb, g, and 
Po- As for volume control, the sedimenting velocity u is 
not a control parameter but determined from demanding 
a force free capsule. 

With the additional pressure parameter pq, we find 
more sedimenting shapes to be stable (with different vol¬ 
umes), in particular, more classes of buckled shapes. 
We tried to find solutions at Eb = 0.01 for pressures 
Po S [—0.95, 0.95] and g G [0,1]. As a consequence of the 
quadratic nature of the elastic law, see Eq. (181, there is 
no solution for a resting capsule {g = 0) with po > 1 [37] . 
For a cleaner comparison of different parameter sets (in 
particular for comparison of differently elongated shapes) 
it is helpful to not use po but the pressure at the center 
of mass, Peff = Po — {zP)g, as an actual control parameter 
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because we use a coordinate system that fixes z = 0 at 
the lower apex of the capsule. 

The possibility to explore a range of volumes allows 
for a greater variety of shapes. In Fig. EH we show 
the pressure-volume relation for stationary axisymmet- 
ric shapes with a velocity close to u = 0.07 (which de¬ 
termines for each shape a certain gravity g). Solutions 
at lower velocities are typically easier to find, but very 
low velocities are rather atypical because one ultimately 
sees the variety of static solutions m- Each type of sta¬ 
ble shape gives one pressure-volume branch in Fig. [m 
Apart from pseudospherical, pear-shaped and the pair 
of weakly and strongly buckled solutions, we find addi¬ 
tional branches, which correspond to higher-order buck¬ 
led shapes with three, four, or five indentations. 

As for fixed volume, the main contribution to the shape 
is due to static pressure (differences). For the pseudo¬ 
spherical branch this can be seen by comparing the vol¬ 


ume (in reduced units) with the known pressure-volume 
relation for a spherical capsule in the static {g = 0) case 

I4phere(p) = ^ ^ (l " ^ (El) 

which fits our data rather well. 

The higher-order buckled shapes occur at small vol¬ 
umes or small effective pressure Peff- We expect to find 
corresponding higher-order buckled shapes also under 
volume control if we prescribe a corresponding small vol¬ 
ume Vb ~ 1 — 2. But also for higher prescribed volumes 
such as Vq = 47 t/3 ~ 4.19, these shapes could become 
stable at very high g: Higher volumes require a larger 
pressure Peff and, thus, also higher deformation forces to 
induce buckling. Such higher deformation forces can be 
achieved at high g. 
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